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Abstract 


Two-dimensional magnetohydrodynamic turbulence is explored by means of 
numerical simulation. Previous analytical theory, based on non-dissipative 
constants of the motion in a truncated Fourier representation, is verified 
by following the evolution of highly non-equilibrium initial conditions 
numerically. Dynamo action (conversion of a significant fraction of tur- 
bulent kinetic energy into long-wavelength magnetic field energy) is ob- 
served. It is conjectured that in the presence of dissipation and external 
forcing, a dual cascade will be observed for zero-helicity situations. 
Energy will cascade to higher wave numbers simultaneously with a cascade 

of mean square vector potential to lower wave numbers, leading to an 

- 1/3 

omni-directional magnetic energy spectrum which varies as k at lower 
wave numbers, simultaneously with a buildup of magnetic excitation at 
the lowest wave number of the system. Equipartition of kinetic and 
magnetic energies is expected at the highest wave numbers in the system. 


2 



1 . 


Introduction 


In a variety of diverse situations, a plasma is known to develop 
magnetic fields on macroscopic spatial scales, with the supporting electric 
current distributions lying interior to the plasma itself. Several 
examples will be familiar: (1) "Magnetic dynamo" generation of the earth’s 
magnetic field; (2) protuberances during solar flares whose geometry 
strongly indicates the motion of magnetic tubes of force; (3) appearance 
of megagauss magnetic fields when solid targets are irradiated by lasers; 

(4) "magnetic island" formation in the poloidal magnetic field of a tokamak 
discharge (though much of this behavior at present is inferred very in- 
directly); and (5) consolidation of electric current filaments generated 
from the Weibel instability in relativistic electron beam simulations. 

Some of these plasma phenomena seem to have little in common, except 
for the decision on the part of the plasma in each case to dump a certain 
fraction of its kinetic energy into the creation of magnetic fields of 
macroscopic dimensions. The processes by which this happens are generally 
poorly understood. It is of interest to search for relatively universal 
mechanisms for generating long-wavelength, self-consistent magnetic 
fields which do not rely upon solving detailed equations of motion for 
highly specialized cases. This is particularly true since many, if not 
all, of the above plasma situations may be highly turbulent. 

A precedent for the conversion of microscopic turbulent energy into 
energy associated with wavelengths comparable in size to the dimensions 
of the system exists in the theory of two-dimensional Navier-Stokes fluids, 
in the form of inverse energy cascades . The cascade of energy to ever- 
shorter wavelengths in three-dimensional Navier-Stokes turbulence is now 
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widely known, as is the inertial-range k kinetic energy spectrum which 
results. Less widely known, but one of the major developments of the 
last ten years in fluid turbulence theory, is the dual cascade which may 
result when boundary and/or initial conditions constrain high Reynolds 
number turbulence to be essentially two-dimensional. (Atmospheric cir- 
culation patterns at scales large comparable with the thickness of the 
earth’s atmosphere are an example.) In this case, evidence is accumulating 
that two quantities are typically cascaded away from any energy source, in 
opposite directions in wave number space. Enstrophy (mean square vorticity) 
has a tendency to flow to higher wave numbers, while energy flows to lower 
wave numbers, thus increasing the degree of long-range organization of the 
flow. Kolmogoroff dimensional arguments applied to the upper and lower 
wave number ranges for steady states give characteristic energy spectra of 
k and k"~ , respectively, together with a continual buildup of energy 

at the longest allowed wavelengths. For analytical, numerical, and expe- 
rimental discussion of these results, see, e.g., Kraichnan (1967); Leith 
(1968); Lilly (1969, 1971); Wiin-Nielsen (1967); Pouquet et al. (1975a); 
and Seyler et al. (1975). 

To the extent that the plasmas in the examples treated in the first 
paragraph of this paper are representable by incompressible magnetohydro- 
dynamic equations, MHD provides a starting point for the dynamo problems 
described in that paragraph. If in some turbulent situations, magnetic 
quantities can be shown to cascade naturally to longer wavelengths, these 
processes will be prime candidates for the development of macroscopic 
magnetic fields. The mathematical structure employed in the Navier-Stokes 
case is readily modified to incompressible MHD. 
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The perception of this possibility appears to have originated with 
the astrophysics group at the Observatoire de Nice, who have concentrated 
on the incompressible three-dimensional case (Frisch et al.» 1975; 

Pouquet et al., 1975b; Pouquet and Patterson, 1976). The present calcu- 
lation being reported here concerns the two-dimensional case, and extends 
an earlier publication by us (Fyfe and Montgomery, 1975). We have chosen 
to investigate the two-dimensional case for three reasons: (1) The two- 
and three-dimensional cases are as different from each other in MHD as 
in Navier-Stokes theory; (2) suggestions have begun to appear that in 
strongly-magnetized situations, MHD turbulence becomes two-dimensional 
with increasing magnetic field strength (Kit and Tsinober, 1971; Schumann, 
1976); (3) an inverse magnetic cascade is possible in two dimensions 
within the framework of a purely isotropic turbulence theory, whereas 
the three-dimensional case requires the presence of magnetic helicity 
and thus a departure from isotropy. It is not known to what extent the 
field variables in the five examples enumerated in the first paragraph 
are satisfactorily representable by situations in which the fields are 
independent of one spatial coordinate, though (1) and (5) surely are, 
and probably others. Finally, external d.c. magnetic fields are easier 
to include in two dimensions than in three; if they are normal to the 
plane of variation of the MHD quantities, they do not affect the mathe- 
matics at all. 

The plan of this manuscript is as follows. Section 2 summarizes 
the results of the inviscid, infinite-conductivity predictions (Fyfe and 
Montgomery, 1975). Section 3, which contains the bulk of the new results 
in this paper, is devoted to the demonstration, by numerical solution of 
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the equations of motion, of the correctness of the statistical-mechanical 
predictions of Section 2. Section 4 discusses cascade processes, which 
necessarily involves the introduction of forcing terms and dissipation 
(finite transport coefficients). Magnetic energy spectra are predicted. 
Numerical simulation to verify these predictions, a considerably more 
ambitious undertaking, is deferred to a third paper in this series. 
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2. Predictions of the non-dissipative theory; relation to cascades 
The first step in discovering a turbulent cascade process can be 
thought of as identification of the rugged constants of the non-dissipative 
motion: these, apparently, are the quantities which can be cascaded. What 
are desired are the integral invariants which result when the dissipative 
terms are dropped from the governing dynamical equations. By "rugged" 
constants of the motion, we mean the sub-class of these invariants which 
are still invariant if both the field variables and the dynamical equations 
suffer truncation in the form of Fourier series representations with large 
but finite numbers of terms. In practice, these are all quadratic inte- 
grals, and many constants of the motion in the untruncated representations 
(e.g., the pointwise vorticity in the two-dimensional Navier-Stokes case) 
are not rugged in this sense. 

Constants of the motion once identified can be used to construct 
thermal equilibrium distributions by the elementary Gibbsian methods of 
statistical mechanics (Kraichnan, 1975), and one can find canonical 
ensembles which have exponential dependences on those constants which 
are "rugged" in the above sense. Multiplying each constant of the motion 
is its associated reciprocal "temperature," which enters as a Lagrange 
multiplier when finding the most probable state. Some of these temperatures 
may be negative, subject only to the overall requirement that the canonical 
distributions shall be normalizable. From these canonical distributions, 
such expectation values as mean energies per Fourier mode are calculable 
in terms of the temperatures by elementary integrations. The temperatures 
themselves are determined by requiring the expectation values of the 
rugged invariants to match prescribed initial values for the system at 
hand. 


7 



Once the rugged constants of the motion are established, the re- 
introduction of dissipation in the system, generally effective only at 
the shorter wavelengths, is necessary for a cascade process and more 
generally, for a truly realistic description of the system. If the 
cascade process is to be seen in its classical, steady-state, Kolmogoroff 
form, a source of the cascaded quantities, or "forcing term," is required 
to balance the dissipation. It is assumed that the forcing terms can 
supply the cascaded quantities at some prescribed rate. Dimensional 
analysis (or more detailed dynamical arguments) lead to wave number 
dependences for the spectra of the cascaded quantities. 

The program described above requires first of all a knowledge that 
all the rugged constants of the motion have been identified. This is in 
some ways the most difficult part of the theory. There is apparently no 
other way to verify that the constants have all been identified than to 
solve numerically the dynamical equations in the truncated Fourier re- 
presentation. The carrying out of this program is the principal content 
of this paper (see Section 3). In particular, thermal equilibrium modal 
energy spectra depend strongly upon which constants of the motion are used 
in the construction of the canonical ensemble. It is of interest, then, 
to follow the evolution of the dynamical equations from highly non-thermal 
initial conditions to verify the spectral predictions for the state to 
which the truncated dynamical equations relax. In Section A, we shall 
return to the completion of the Kolmogoroff program for the cascade pre- 
dictions. 

We now summarize the results of the two-dimensional incompressible 
equilibrium theory (Fyfe and Montgomery, 1975). We consider two-dimensional. 
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ideal magneto-hydrodynamics in the case in which the field variables have 

the following geometry in Cartesian coordinates. B » (B^, B y , 0) is the 

magnetic field; v ■ (v , v , 0) is the fluid velocity; j - VxB “ (0, 0, j ) 

■. x y — — 

is the electric current density, u * Vxv ■ (0, 0, to ) is the vorticity 
vector, a “ (0, 0, a ) is the vector potential, and Vxa - B. E « <0, 0, E ) 

2 % % ^ Z 

is the electric field. All quantities are independent of the z coordinate, 
but depend upon x, y, and the time t. A uniform d.c. magnetic field 
B * (0, 0, B.) can be added at any step of the development without changing 
the mathematics. 

All quantities are Fourier-decomposed assuming periodic boundary 
conditions over a large square box; e.g.. 


j - E j (k,t) exp (ik • x) 
z k ~ ■ 


to ■ E to(k,t) exp (ik • x) (1) 

z k - - ' 


and so on. The k » (k , k , 0) are summed over all the values (2‘n , /L)* 

— ^ y 

2 2 

(n , n , 0) where n , n are any integers for which n + n >0, and L is 
x y x y x y 

the linear dimension of the box. 

In the Fourier representation, and in the appropriate dimensionless 
units, the dynamical equations can be written entirely in terms of to(k,t) 
and j(k,t) (Fyfe and Montgomery, 1975): 


3to(k,t) 

9t 

3j(k,t) 

“at 


S M.(r,p) 6 (p + r-k)[to(r)u(p) - j (r) j (p) ) 
2 M 2 ( 5»E> 6 (p + 5-k)[j(r)to(p) - to(r)j (p) ] 


( 2 ) 

(3) 
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where the coupling coefficients and independent of the field 
amplitudes, are 


M,(r»p) - M. (p,r) - (1/2) e «(rxp) 

1 1 ^ Z 


- (? ■ ?) 


( 4 ) 


and 


M 2 (r,p) » -M 2 (p,r) " C 1 / 2 ) ® 2 *^ X E^ ~ 2~2 

P r 


(5) 


Here, e w is a unit vector in the z-direction, and r»p are dummy wave 

numbers which lie in the xy plane and which are summed over all the 

values allowed by the periodic boundary conditions. 6(p + r - k) ** 1 if 

p + r « and is zero otherwise. 

Three invariants have been found which survive the restriction of 

the allowed wave numbers in Eqs. (1), (2), and (3) to lie inside a 

circular annulus between Ik I * k , and jkl ■ k . k . is a minimum 

wave number determined by the box size, equal to 2 tt/L, and k » k , 

max mm 

is a large but finite maximum wave number which, within the framework of 
the non-dissipative theory, remains physically undetermined. These 
three invariants are the total energy £, the cross helicity P, and the 
mean square vector potential A. Referred to unit volume, they are 


€ = *| IT 2 / dxdy (v 2 + B 2 ) ■ y ^ 


|u>(k)| 2 + |j(k)| 


( 6 ) 


• i -‘7 


dxdy v • B 


IE 


u(k) j(-k) 


(7) 
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( 8 ) 



(We have suppressed the time arguments for economy of notation*) 

The canonical distribution, in the phase space defined by the real 
and imaginary parts of the j (k) and w(k) , is 

D « r> exp 2 (-a e-S? -yA) (9) 

wnete n is a normalizing constant, and 3 *,B \y ^ are the energy, 
helicity, and vector potential temperatures, respectively. 

Throughout the balance of this paper, we restrict ourselves to non- 
helical flows, for which <P> * 0; from this It can be shown (Fyfe and 
Montgomery, 1975) that S ■ 0. The case of helical flows is deferred to 
a later publication. Omitting helicity also permits a wholly isotropic 
formulation of the theory. 

For the case of zero helicity, the expectation values calculated 
from Eq. (9) are 


< |B(k)| 2 > - (a + Yk'V 1 

•m <■» 


< 10 ) 


and 


< | v 00 1 2 > - a* 1 . (11) 

•** 

The a and y are determined as solutions of the pair of simultaneous 
equations 
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c 


( 12 ) 


% 


1 

2 


V J(a + Y k-V 1 + a- 1 | 

k ' 1 


and 


| J (V + k 2 a) _1 ■ A, (13) 

k 

where C and A are the initial values of the energy and Integrated, squared, 
vector potential for the system under consideration. 

More convenient than Eqs. (10) and (11) for numerical investigation 
are the equivalent spectra for stream function and vector potential 
given, from Eqs. (10) and (11), by 


and 


l2 < l v < k >l > 21 

< Mk)| Z > “ ^ - (Ctk 2 )' 1 


(14) 


. .2 <|B(k)r> 2 , 

< |a(k)| > - r - (Y + ak Z ) 1 . 


(15) 


a has been shown to be necessarily positive. But Y can either be 

positive or negative. For fixed C and A, y always becomes negative if 

k is taken large enough. As k gets larger and larger, Y + ctk 2 , 
max max min 

approaches zero, and the vector potential becomes more and more concentrated 
in the lowest wave number of the system. This is the first and crucial 
indication of a possible inverse cascade of mean square vector potential. 

It should be noted that the sums in Eqs. (6), (7), (8) are over all 
allowed wave numbers in all four quadrants in k. Not all of these k's 
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correspond to independent Fourier coefficients, since the Fourier coefficients 
in half of the k space determine chose In the other half by the reality 
conditions w(k) * u (-k) and j(k) - J -k). For numerical purposes it Is 
more convenient (and equally acceptable) to carry out those suras only over 
the region k x £ 0, omitting the nodes for which k x - 0 and k y < 0, in 
Eqs. (6), (7), (8), (12) and (13). That is, we may if ve choose sura only 
over tha independent k-modes. This is done in Sec. III. 

•m 
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3. Numerical simulation, and comparison with theory 


The d ^namical equations are solved by a spectral method closely 
related to the method of Orszag (1971) and Patterson and Orszag (1971) 
(see also: Salu and Knorr, 1975 and Sevier et al., 1975). Rather than 

solving Eqs. (2) and (3), which are more convenient for theoretical 
purposes, it is preferable numerically to solve the equivalent pair of 
equations for the vorticity and vector potential. 




a -vVuj + B»V1 
ct z z 

(16) 

3a 


ar ■ -r 7 s 

(17) 


with the various quantities appropriately Fourier-decomposed. All other 
field quantities are expressible in terms of the solutions tc (16) and 
(17). It is of the essence of the spectral method that the products on 
the right hand sides of Eqs. (16) and (17) are evaluated in position 
space rather than as convolution sums, with a fast Fourier transform 
back to k-space before the Fourier coefficients are advanced at the next 
time step (see Patterson and Orszag, 1971 or Salu and Knorr, 1975). 

Initial conditions are specified by giving the t = 0 values of the 
Fourier coefficients. This also facilitates numerical solution of 
Eqs. '12) and (13) for a and y. As in our earlier work (Seyler et al. , 
1975), we choose highly non-thermal initial conditions, by setting all 
Fourier coefficients initially zero outside a ring in k-space. All 
excitations then move to both larger and smaller values of |k| according 
to the solutions of (16) and (17). Eventually, a state is reached in 


which the amplitudes, though still fluctuating with time, show no systematic 

drift toward either larger or smaller values. We then time-average the 

squares of the Fourier coefficients over a long enough interval that the 

average becomes essentially independent of time. The upper limits on 

the time intervals over which we may average are indicated by the time 

over which our three constants of the motion (C» A, and P = 0) drift 

away from their initial values. These time averages are then compared 

with the theoretical predictions of Eqs. (14) and (15). 

Several runs were carried out, and we illustrate the results by a 

description of the two runs to which the greatest amounts of computer 

time were devoted, and upon which the most diagnostics were performed. 

One of these, called run Fl hereafter, corresponds to y < 0; the other, 

called run J2 hereafter, corresponds to y > 0. Table 1 summarizes the 

important parsm rvrs for these two runs. 

Fig. 1 consists of three panels. Fig. la shows the initial values 
2 2 

of |a(k)| vs. k for run Fl. Fig. lb shows the time averages of the 
2 

|a(k)j between times 17.5 and 26.25 for run Fl. Fig. lc shows the time 

2 

average of the ja(k)| between times 35.0 and 52.5, at the end of run Fl. 

The solid curve in Fig. lc is the theoretical curve, Eq. (15). Different 

2 

values of k which correspond to the same value of k have been averaged 
over in all spectral plots. 

Fig. 2 consists also of three panels, which are the corresponding 

2 

stream function spectra t^j(k)| for run Fl. Fig. 2a shows the initial 

2 2 
j^(k)| , Fig. 2b shows the time averaged |i^(k)[ between times 17.5 and 

*"W “* 

2 

26.25, and Fig, 2c shows the time averaged |^(k)| between times 35.0 
and 52.5, at the end of the run. The solid curve in Fig. 2c is the 
theoretical curve, Eq. (14). 
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Figs. 3 and 4, consisting of three panels each, show the corresponding 

quantities for run J2. Fig. 3a is the initial data for the |a(k)| 2 , 

2 

Fig. 3b is a time average of the ja(k)|“ between times 24.0 and 31.0, 

2 

and Fig. 3c is the time average of j a ( k) 1 between times 34.5 and 51.0, 

at the end of run J2. The solid line in Fig. 3c is the theoretical 

curve, Eq. (15). The three panels of Fig. 4 are the corresponding 

2 2 

stream function plots (j^(k)| vs. k ) for run J2. 

Figs. 1 through 4 are considered to be verification of the spectral 
predictions of the Eqs. (14) and (15), the most important aspects of the 
non-dissipative equilibrium theory. 

Dynamo action (i.e., the conversion of kinetic energy into magnetic 
energy) is shown in Fig. 5. Fig. 5 is a plot of the ratios of total 
magnetic energy to total kinetic energy versus time, for both runs. The 
theoretical predictions are indicated by the arrows on the right hand 
side of the graph. 

The time history of representative modal energies, averaged over 
2 

inodes for a fixed k , but differing k, are shown in Fig. 6 for run FI. 

1 2 2 

Fig. 6a shows } B (k) | vs. time for k * 1, 2, 8, and 26. Fig. 6b shows 

2 2 

| v (k) j vs. time fcr the same values of k . It is characteristic of the 

•V *v 

magnetohydrodynamic problem (and is the principal observed difference 
from the two-dimensional Navier-Stokes case) that the temporal variations 
of the | B (k) 1 *" for low values of k are significantly slower than the 
temporal variations of the |v(k)| , or of the |s(k)| for higher values 
of k. This confirms qualitatively an expectation raised by noting that, 
in Eqs. (2) to (5) , 
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Hi ( r * p) 

for k « p. Thus the flow of magnetic energy to the lowest order modes 

is expected to be very slow compared to the flow of kinetic energy among 

the same modes. This can be seen in Fig. 6. 

Fig. 7 shows contour plots of constant a at four different times 

z 

in the evolution of run FI. Fig. 8 shows similar contour plots for 
constant stream function y . Fig. 9 shows contour plots of constant a 

z 

at four different times during the evolution of run J2, and Fig. 10 

shows the evolution of the stream function ^ for run J2 in the same 

fashion. Particularly in Fig. 7, it will be seen that the migration of 

magnetic energy to longer wavelengths is accompanied by the growth and 

consolidation of the contours of constant a (which, in this geometry, 

z 

are also magnetic field lines). The measure of the total area occupied 
by any range of values remains, however, approximately constant. 

The result has been often quoted, but apparently not proved, that 
ideal magnetohydrodynamics permits no changes in the topologies of 
closed magnetic field lines in two-dimensional situations in magneto- 
fluids. We do not know if the result is true, but it should be kept in 
mind that in any case, the result certainly does not apply to two- 
dimensional MHD in the truncated Fourier representation. If a proof 
could be constructed, it would undoubtedly rely on the invariance of the 
two-dimensional analogue of the Alfven flux invariant, 

2 

I » I e^ ■ J" di x B, 

1 
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where d£ is a vector element of length along a curve moving with the 

fluid between end points 1 and 2. While I is invariant for the original 

ideal MHD equations before Fourier decomposition, it does not survive 

the restriction of Eqs. (2) and (3) to large but finite numbers of 

terms. I is, like the pointwise vorticity in the ideal Navier-Stokes 

case (and like a in the ideal two-dimensional MHD case) not a "rugged" 

invariant in the sense of Sec. II. Therefore, no theorems based on its 

invariance can apply to the present system. This leaves open, as the 

question remains open in the Navier-Stokes case, the question of the 

relation of the truncated system (Eqs. (2) and (3)) to the original 

system of ideal MHD equations. This is an important unanswered question, 

though its importance is perhaps diminished by the addition of finite 

viscosity and resistivity (Sec. IV), which likewise destroy the invariance 

of I. No mathematical investigation of the ideal incompressible MHD 

equations carried out so far begins to approach in rigor what would be 

necessary for a resolution of this question. For example, possible 

singularities which may develop in the ideal MHD equations in the incompressible 

case have not been addressed; investigations are still stalled for the 

case of ideal Navier-Stokes fluids, which are a special case of ideal 

MHD. For this reason it seems necessary to leave this question open for 

the time being. 

That "magnetic islands" do indeed seem to merge in the present 
truncated representation can be readily seen by looking at sequences of 
typical closely (temporally) spaced a^-contour plots. Some of these are 
shown in Fig. 11. During the sequence shown in Fig. 11, £ and A vary by 
only .05% and .01%, respectively, so it is not likely that numerical 
inaccuracies constitute a "dissipation" which is sufficient to overcome 
any real or conjectured limitations ideal MHD may place on the merging. 
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4 . Addition of dissipation and forcing ; spectral predictions 

The presence of finite dissipation and forcing modifies Eqs. (2) 
and (3) to read 


t>u>(k,t) 

It 


^ M 1 (r,p)5 (p+r-k) [o>(r)u)(p)-j (r)j(p)] 


- vk 2 U)(k,t) + f(k,t) 


aj(k.t) 

^ * 2J M 2 (£»P)<5 Cp-*-£-k) [ j (r)oj(p)Ho(r)j(p)) 


(18) 


- yk 2 j(k,t) + g(k,t) . 


(19) 


In Eqs. (18) and (19), v is a dimensionless viscosity and y is a dimen- 
sionless resistivity. f(k,t) and g(k,t) are prescribed functions of the 
time which model external forces (mechanical and/or magnetic) which may 
act on the system. Numerical investigation of the system (18) and (19) 
is planned. At present, we are not interested in detailed solutions of 
(18) and (19), but in generalizing the Kolmogoroff program of dimensional 
analysis to make predictions about possible cascades of G and A which 
Eqs. (18) and (19) may imply, when v and y are very small but non-zero. 

Cascade processes have been considered from two different points of 
view: the original, relatively simple Kolmogoroff (1941) dimensional 

arguments, and much more sophisticated dynamical arguments which lead, 
e.g., to logarithmic corrections (cf. Kraichnan, 1973) and which depend 
upon detailed properties of the Navier-Stokes equations. At the outset, 
it seems desirable to stay with the simpler point of view, since little 
is as yet available in the way of accurate information on the detailed 
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properties of the l!HD turbulent solutions, and since it is as yet unclear 
what the increased sophistication has added to the Kolmogoroff result. 

A natural generalization of Kolmogoroff 's program seems to be to 
construct cascadable quantities per unit mass by dividing extensive 
rugged invariants by the total mass of the system. There are two of 
these for this situation, the energy per unit mass and the integrated 
vector potential per unit mass. We abandon the dimensionless variables 
at this point, which can obscure the dimensional analysis, and revert to 
c.g.s. laboratory variables. The quantities to be cascaded are therefore 
(p * mass density) 


n = 


a = 


Jib: 

f 

/ e 
/ 


pv^ + B^/8irj dxdy 


p dxdy 


a^ dxdy 


p dxdy 


( 20 ) 


( 21 ) 


2 -2 

H is the total energy per unit mass and has dimensionality L T , while 

A is the integrated vector potential per unit mass and has dimensions 
A 2 

L T . It is assumed that external sources are capable of supplying II 
and A to the medium at the rates £ per unit time and h per unit time 
respectively. 

First, a decision must be made in any dual-cascade situation as to 
which quantity will cascade in which direction. We proceed by analogy 
with the Navier-Stokes case, and assume that the quantity which in the 
the non-dissipative equilibrium theory can be sharply-peaked at small k 
is the one which can be cascaded to longer wavelength: namely A. It is 

also reasonable that the cascaded quantity which contains the highest 
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power of k in its defining expression should be cascaded to higher k, 
since it is the easier quantity to dissipate. The spectrum (omni- 
directional spectrum) A^> related to A by 

CO 

/ A k dk = * 

0 

5 -2 

has dimensionality L T . If is assumed to depend only upon n and k 

OL 8 

in the combination n k , simple arithmetic shows that 


a - 2/3 
8 - -7/3 . 


( 22 ) 


The essential content of the second of Iqs. (22) is that the omni- 

directional magnetic energy spectrum , 2“ k< ! B (k) | >, should vary as 
-1/3 

k . Of course, eventual pile-up of the magnetic energy at the longest 
wavelength will excite the lowest modes even above this level, as in the 
Navier-Stokes case. 

Eqs. (22) apply below the input wave numbers characterizing the 
sources f(k,t) and g(k,t). Above these input wave numbers, similar 

-5/3 

dimensional arguments predict an omni-directional energy spectrum •"k 

Kraichnan (1965), (see also Pouquet et al. , 1975b) has proposed in three 

dimensions an "Alfven effect” which inhibits energy transfer at the 

higher wave numbers and would lead to an omni-directional energy spectrum 
-3/2 -5/3 

-k instead of ~k . Present numerical experiments may find it 
difficult to distinguish between these two possibilities. In any case, 
it is our opinion that the question of most interest is the possibility 
of magnetic excitations at the longest wave numbers rather than the 
precise value of the exponent at large wave numbers. It is this effect 
which leads to the dynamo action described here. 
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5. Discussion 


It Is considered that two basic contributions are contained in the 
preceding pages: (1) Verification to a considerable extent of some 

previous predictions about equilibrium statistics of a truncated Fourier 
representation of the two dimensional MHD equations; (2) A conjecture 
concerning the cascade of vector potential to longer wavelengths for the 
same system In the presence of dissipation and external forcing. 

It is important to Isolate, also, two or three aspects of the 
problem that have not been resolved and which require further study. 

Most important of these (which parallels long-standing problems in 
Navier-Stokes turbulence theory) is the relation of the original ideal 
MHD equations to their representation in terms of a large but finite 
Fourier representation. Experience with other systems of ideal fluid 
equations for which exact information is available suggests the possibility 
that singularities may form after finite times for most sets of initial 
conditions. If this is the case, these singularities will not be visible 
in any finite Fourier representation. It is also possible that while 
such singularities might form for the ideal MHD equations, the addition 
of small but finite dissipation might domesticate them (e.g.. Grad, 

1975). In this case also, the relation of the truncated Fourier repre- 
sentation might obscure the process to some degree. However, all these 
objections can be raised in simpler form in connection with Navier- 
Stokes turbulence, and the fact that they are unanswered there has not 
altered the fact that virtually everything we know about Navier-Stokes 
turbulence has come about through truncated Fourier representations. To 
those who are to some degree familiar with the current status of Navier- 
Stokes turbulence theory, the above limitations may seem less serious 
than to those who are not. 
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Table 1 



Run FI 

Run J2 

Independent k values for which 
u(k) ^ 0, initially 

(3,1) , (3,2) , (2,3) 
(3,-1) , (3,-2) , (2,-3) 
(1,3) (1,-3) 

same as 

FI 

Independent & values for which 
j (k) j 0, initially 

(2,2), (3,0) 
(0,3), (2,-2) 

same as 

FI 

Initial C 

4.002 

3.962 

Initial A 

0.125 

0.0484 

Initial P 

zero 

zero 

a (from initial data) 
a (from final data) 

183.8 

169.5 

179.6 

169.0 

Y (from initial data) 

Y (from final data) 

-157.4 

-142.0 

92.96 

135.9 

Duration of run (t * t ) 

max 

52.5 

52.5 

Size of time step (it) 

(128)~ X 

(128)' 1 

Final *£ 

4.329 

4.197 

Final A 

0.126 

0.0490 

Final P 

0.065 

-0.0112 

Percent change in C (at end of run) 

8.17 

5.93 

Percent change in A (at end of run) 

0.68 

1.23 

k . - 2 tt/L 

min 

1 

1 

k 

max 

16 

16 

(Magnetic energy/kinetic energy) Ratio 
Initial 

Average over last 2240 time steps 
Theoretical (from initial data) 

0.361 

1.047 

1.055 

0.1159 

.9857 

.9873 


REPRODUCIBILITY OF THE 

’ * A ' n< ’’> TR POOR 
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Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Fig. 6 


Fig, / 


Figure Captions 

2 2 

Vector potential spectra ja(k)| vs. k for run FI: (a) initial 
conditions; (b) average over tines t - 17.5 to 26.5; (c) aver- 
age over times 35.0 to 52.5. Theoretical curve Is solid line 

2 

in Fig. lc. For clarity, only every other value of k is 

2 2 
plotted for 50<k <100 and every fourth value for k >100; the 

same convention is adopted also in Figs. 2 , 3, and 4. 

2 

Stream function spectra |^(k)| for run FI: (a) initial conditions; 

(b) average over times 17.5 to 26.5; (c) average over times 35.0 

to 52.5. Theoretical curve is solid line in Fig. 2c. 

2 2 

Vector potential spectra |a(k)| vs. k for run J2: (a) initial 

•» 

conditions; (b) average over times 24.0 to 31.0, (c) average over 
times 34.5 to 51.0. Theoretical curve is solid line in Fig. 3c. 
Stream function spectra jo(k)J vs. k for run J2: (a) initial 
conditions; (b) average over times 24.0 to 31.0; (c) average 
over times 34.5 to 51.0. Theoretical curve is solid line in 
Fig. 4c. 

Computed ratio of total magnetic energy to total kinetic energy 
versus time, for runs Fl end J2. Theoretical prediction is 
Indicated at right of graph. 

2 

Time history of ene *ies for fixed k ■ 1, 2, 8, and 26. 

2 2 
(a) ]B(k)j , magnetic energy, and (b)|v(k)| , kinetic energy. 

Both are for run Fl. 

Contour plots * const at four different times in the 
evolution of run Fl. 
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Fig. 

Fig. 

Fig. 

Fig. 


S Contour plots y * const for stream function at four different 
tim*s for run FI. 

9 Contour plots «. * const at four different times during the 

evolution of run J2. 

10 Contone plots v ■ const at four different times during the 
evolution of run J2. 

11 Four contour plots a z ■ const for run FI at closely spaced times 
showing growth and consolidation of magnetic contours, 
characteristic of y < 0. 
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Figures la, 



REPRODUCIBILITY of the 
ORIGINAL PAGE SB POOR 


Figures 







Figures 3a 









Figures 4a , 









MAGNETIC ENERGY/KINETIC ENERGY 
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REPRODUCIBILITY OF THF 
OTUOJNAL PACE IS POOP 






figure 8 



















